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Abstract 


An  investigation  is  made  as  to  whether  or  not  pulse 
shapes  within  a  single  planar  Ge(Li)  detector  can  be 
used  to  measure  gamma  ray  polarization,  Compton 
scattering  polarization  effects  are  studied  experimentally 
and  theoretically  using  both  a  Triple  Constant  Fraction 
Pulse  Shape  Discriminator  and  a  differentiating  pulse 
shape  analyser.  As  well,  a  short  theoretical  study  is 
made  of  photoelectric  absorption  polarization  effects. 

The  results  of  these  studies  show  that  pulse  shape 
analysis  in  a  Ge(Li)  detector  is  not  a  practical  method 
of  measuring  gamma  ray  polarization. 
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Chapter  I 


Introduction 

During  the  past  decade  considerable  progress  has 
been  made  in  the  development  of  nuclear  semiconductor 
detectors.  The  introduction  of  large  volume  lithium 
drifted  germanium,  or  Ge(Li),  detectors  has  provided 
experimenters  with  an  instrument  for  gamma  ray  detection 
which  combines  both  excellent  energy  resolution  with 
reasonable  overall  detection  efficiency.  The  use  of 
this  device  has  therefore  led  to  advances  in  many  areas 
of  nuclear  gamma  ray  studies,  including  that  branch  which 
deals  with  gamma  polarization  measurement.  As  a  result, 
in  recent  years  a  number  of  gamma  polarimeters,  all 
capable  of  good  energy  resolution,  have  been  developed. 

These  polarimeters,  which  make  use  of  the  scattering 
asymmetry  inherent  in  the  compton  scattering  process, 
have  included  both  conventional  two-  and  three-crystal 
compton  polarimeters  and  single  crystal  planar  detectors. 
The  planar  detectors,  in  particular,  appear  to  be  quite 
attractive,  due  to  their  simplicity  of  experimental 
setup.  Not  only  are  they  simpler  physically  to  arrange 
and  manipulate,  but  also  they  are  less  complicated 
electronically  to  work  with,  since  they  do  not  require 
much  of  the  excess  electronic  equipment  necessary  to 
handle  a  multidetector  coincidence  system.  However, 
recent  studies  have  shown  that,  although  the  planar 
detector  has  a  higher  overall  detection  efficiency  than 
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multicrystal  systems,  its  relative  polarization 
sensitivity  is  quite  small.  Hence  it  would  be  of 
considerable  benefit  to  design  a  single  crystal  gamma 
detector  capable  of  high  polarization  sensitivity. 

This  thesis  investigates  the  question  of  whether 
or  not  pulse  shapes  within  a  single  detector  can  be 
used  to  measure  the  polarization  asymmetry  produced  by 
compton  scattering.  Two  different  methods  of  measuring 
these  shapes  are  explained,  and  their  individual  merits 
and  deficiencies  discussed.  As  well,  a  short  appendix 
is  included  which  studies  the  polarization  effect 
inherent  in  the  photoelectric  absorption  of  gamma  photons. 


Chapter  II 


Theory 

2.1  Compton  Scattering;  and  Polarization  Measurement 

It  can  be  shown  that  the  differential  collision 
cross  section  for  the  scattering  of  incident  plane- 
polarized  radiation  is  given  by 

is-  __  (vi  +  Y_'  -  Z  +  H  co%1  x 

AR-  H  \VS  I  \V  K 

where  X  is  the  angle  between  the  electric  vectors 
of  the  incident  and  scattered  radiation,  hVo  is  the 
energy  of  the  incident  photon,  Ky7  the  energy  of  the 
scattered  photon,  and  r0  the  classical  electron  radius. 
The  angles  and  directions  encountered  in  the  compton 
scattering  process  are  visualized  in  fig.  1,  in  which: 
OA  ®  direction  of  scattered  photon, 

OAB  =  scattering  plane, 

■©•  =  scattering  angle, 

&o  =  electric  vector  of  incident  photon, 

^  =  angle  between  60  and  OA, 

ODB  «  plane  of  polarization  of  incident  radiation, 

4>  =  angle  between  plane  of  incident  polarization  and 
scattering  plane, 

=  electric  vector  of  scattered  photon, 

€,j  =  projection  of  £7  on  OAC  plane, 

£l  =  projection  of  £*  perpendicular  to  OAC  plane, 
p  =  angle  between  £/  and  OAC  plane,  and 
X  =  angle  between  £  and  E0. 
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Figure  1  Compton  Scattering  Terminology 
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The  energy  of  the  scattered  photon  is  given  by: 

i  ?  i 

n  y  -  m0 c 


I  -co^e-t  ( \/oc) 


where 


oC  z  KV0 


Thus 


Thus  .  „ / 


) 


Vo  I  -^cxCi-cos^) 

By  geometrical  arguments  it  can  be  shown  that 


cosX  *  cos^  sin^ 


and 


Therefore,  summing  over  all  possible  directions  of  final 
polarization, 


Hence  it  is  apparent  that  plane-polarized  radiation 
will  tend  to  be  scattered  in  a  plane  perpendicular  to 
the  plane  of  polarization  of  the  incident  radiation 


(<j>  =  90°). 


This  polarization  effect  is  used  as  the  basis  for 
many  polarization  measuring  devises,  including  the  two- 
and  three-crystal  compton  polarimeters  and  the  single 
crystal  planar  detector.  The  standard  experimental 
procedure  involves  the  measurement  of  a  scattering 
asymmetry  A,  defined  by  A  =  QP,  where  P  is  the  linear 
polarization  of  the  incident  photon  beam  (-1  -  P  £  +1), 
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and  Q  is  the  polarization  sensitivity  of  the  detector 


(0  -  Q  £  +1).  "A"  is  measured  by  comparing  the  counting 


rates  of  the  polarimeter  when  oriented  in  the  U)  *  0° 


direction  and  the  CJ  -  90°  direction  (see  fig.  2). 
Thus 


where  N((J)  =  the  number  of  counts  occuring  in  the  detector 
when  it  is  oriented  at  the  angle  UJ. 

There  are  two  ways  this  asymmetry  in  scattering 
can  be  detected  when  only  using  a  single  detector.  To 
begin  with,  the  detector  crystal  can  be  of  some  shape 
which  preferentially  absorbes  photons  scattered  in  a 
particular  direction,  as  is  the  case  for  a  planar  detector. 
Unfortunately,  the  polarization  sensitivity  of  this 
experimental  setup  is  quite  poor,  even  though  its  overall 
detection  efficiency  is  high.  Hence  it  is  normally  used 
only  in  cases  where  the  overall  accuracy  of  a  polarization 
measurement  is  limited  by  its  counting  statistics. 

Secondly,  though,  the  scattering  asymmetry  could  be 
measured  within  a  single  crystal  by  considering  the 
shape  of  the  electronic  pulses  triggered  within  it.  To 
understand  how  this  could  be  done,  one  must  first  under¬ 
stand  how  electronic  signals  are  formed  in  Ge(Li) 
detectors,  and  as  well  how  they  can  or  could  be  measured. 

2.2  Pulse  Shapes  in  a  p-i-n  Detector 


The  simplest  way  to  determine  the  pulse  shape 
produced  by  a  planar  Ge(Li)  detector  is  to  first  consider 
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Figure  2  Experimental  Asymmetry  in  Gamma  Ray  Polarimeters 
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the  shape  produced  when  only  a  single  electron-hole 
(e-h)  pair  is  created  within  the  detector  crystal. 

This  situation  is  depicted  in  fig.  3a,  in  which  d  is 
the  detector  thickness,  V  the  voltage  drop  across  it, 
and  Z©  the  distance  from  the  creation  point  of  the  e-h 
pair  to  the  crystal's  p  layer.  If  jUz  and  are  the 
electron  and  hole  mobilities,  then  the  currents  in  the 
external  circuit  due  to  the  electron  and  hole  are 
ie  =  (V/d2)  q  yUe 

and  ih  =  (V/d2)  q 

where  of  course  q  is  just  the  elementary  charge  of  an 
electron.  The  time  taken  for  the  electron  and  hole  to 
reach  the  outer  layers  of  the  detector  is 
Te  =  {(d  -  (d/V) 

and  Th  =  {Z./yU^  (d/V)  . 

Since  the  final  output  signal  is  just  the  sum  of  the 
integrated  electron  and  hole  current  signals,  the  total 
signal  out  of  the  detector  is  as  shown  in  fig.  3b.  Note 
that  the  total  charge  flow  in  the  external  circuit  is 

^total  ~  ^e^e  +  ^h^h 

=  q  • 

Most  germanium  detectors,  and  in  particular  Ge(Li) 
detectors  are  operated  at  a  temperature  of  about  77°K. 

At  this  temperature  the  electron  mobility  in  germanium 
is  approximately  equal  to  the  hole  mobility,  and  thus 
ie  =  ih*  Assuming,  therefore,  that  the  electric  field 
inside  the  planar  Ge(Li)  is  in  fact  a  constant  throughout 


■ 

& 
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the  crystal,  and  that  electronic  noise  in  the  system  is 
minimal,  the  total  output  signal  for  a  single  e-h  pair 
will  appear  as  in  fig.  3c.  It  should  be  noted  here 
that  the  shape  of  the  signal  is  dependent  upon  the  location 
of  creation  of  the  e-h  pair. 

Of  course  all  of  the  above  discussion  is  quite 
simplified  when  compared  to  what  actually  happens  in  the 
lab.  During  a  gamma-ray  detection  experiment  a  single 
photon  may  be  photo-absorbed,  compton  scattered,  or 
absorbed  by  pair-production  within  the  Ge(Li)  crystal. 

At  each  point  in  the  detector  where  the  gamma  loses  some 
energy,  a  large  number  of  e-h  pairs  are  produced.  Further¬ 
more,  the  ’’point"  at  which  they  are  produced  has  finite 
dimensions,  which  depend  both  on  the  gamma  energy  lost 
there  (see  fig.  4)  and  the  scattering  kinematics  of  the 
ionized  electrons.  Thus  the  real  pulse  shapes  produced 
by  a  planar  Ge(Li)  detector  are  similar  to  those  shown 
in  fig.  5* 
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Figure  4  Electron  Range  as  a  Function  of  Energy 
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Figure  5  Photographs  of  pulse  shapes  obtained  using  a 
planar  Ge(Li)  detector  and  a  ^°Go  gamma  ray 
source.  The  particular  preamplifier  used  with 
this  detector  provided  a  negative  output  signal. 
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Chapter  I  IT 


Measurement  of  Pulse  Shapes 

3*1  The  Triple  Constant  Fraction  Pulse  Shape  Discriminator 

Theory: 

Once  the  general  appearance  of  the  pulse  shapes  is 
known,  the  question  arises  as  to  whether  or  not  these 
pulse  shapes  can  be  "measured",  and  if  indeed  these 
measurements  are  sensitive  to  the  polarization  effect 
discussed  earlier.  One  device  which  has  been  used  in 
the  past  to  measure  pulse  shape  is  the  triple  constant 

-I 

fraction  pulse  shape  discriminator  (TCF).) 

The  TCF  makes  use  of  two  constant  fraction  timing 
triggers  (CFT)  to  measure  pulse  shape.  Each  CFT  sees  an 
input  pulse  with  a  leading  edge  of  some  shape  g(t),  and 
inverts  and  shrinks  it  in  amplitude  by  a  fraction  f^  (see 
fig.  6).  The  CFT  then  combines  this  smaller  pulse  with 
the  original  pulse,  which  has  been  delayed  by  a  time 
interval  "a",  and  triggers  when  the  resultant  pulse 
crosses  zero.  This  occures  at  a  time  t^ ,  such  that 
g(t1)  =  f1g(t1  +  a)  . 

Note  that  for  a  linear  pulse 


g(t)  »  mt 


Thus  t^  is  independent  of  the  slope  m.  In  a  TCF  one 
CFT  is  set  at  a  fraction  f^  and  used  to  trigger  a  start 
signal  in  a  time-to-analog  converter  (TAC),  while  another 
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Original  Pulse  g(t) 


Figure  6  Triggering  In  A  CFT 
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OFT  is  set  to  a  larger  fraction  f^  and  used  to  trigger 
a  stop  pulse  in  the  TAG.  This  then  defines  a  time  interval 
such  that  &  =  t2  -  t^  where 
g(t^)  =  f1g(t1  +  a) 
and  g(t2)  =  f2g(t2  +  a). 

h  can  therefore  be  used  as  a  rough  measurement  of  the 
pulse  shapes.  For  instance,  note  that  for  a  pulse  with 
a  linear  leading  edge: 

&  =  a  f «  —  f  > 

(i-hXi-A  ) 

=  constant. 

Of  course  the  interpretation  of  &  becomes  much 
more  difficult  when  real  pulse  shapes  of  the  kind  shown 
in  fig.  5  are  encountered.  S  is  highly  dependent  upon 
the  type  and  energy  of  the  interaction  encountered,  the 
position  or  positions  of  the  interaction,  and  in  partic¬ 
ular  the  delay  setting  "a”  of  the  TGF.  Figs.  7  and  8 
show  the  various  i  which  would  be  expected  for  photo¬ 
electric  and  one-compton  events  occuring  at  various 
locations  throughout  the  crystal.  In  these  graphs: 
a  =  the  delay  setting  on  the  TGF, 

=  the  location  of  the  outermost  gamma  interaction, 

(0  refers  to  the  center  of  the  crystal,  1  to  the 
detector  edge), 

Z2  =  the  location  of  the  innermost  gamma  interaction,  and 
=  the  amount  of  the  total  gamma  energy  deposited  at 


location  Z^ 


. 


0.4- 


Photoelectric 

Interaction 


^-o 


1.0 


Figure  7  Delta  as  a  function  of  gamma  interaction 

type  and  location  (a  =  0.5  arbitrary  units) 
Delta  is  measured  in  arbitrary  units. 
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Figure  8  Delta  as  a  function  of  gamma  interaction 

type  and  location  (a  >  2.0  arbitrary  units). 


' 


For  the  purpose  of  these  graphs,  all  pul3e  shapes  were 
considered  to  be  sharp;  that  is  the  range  R0  of  all 
scattered  electrons  was  assumed  to  be  zero. 

Results : 

The  shape  of  a  S  spectrum  and  its  sensitivity  to 
polarization  were  first  investigated  experimentally  using 
a  thin  planar  detector.  This  detector  had  a  radius  of 
1.73  cm  and  a  thickness  of  0.6  cm.  In  the  experiment 
3*4  MeV  protons  were  used  to  excite  ^Fe  nuclei  to  their 
first  excited  state,  from  which  they  decayed  emitting 
photons  of  energy  0.847  MeV.  Angular  correlation  measure¬ 
ments  showed  that  this  radiation  had  a  polarization  at 
90°  of  (0.51*0.03),  while  at  0°  its  polarization  was 
zero.  Shape  spectra  were  accumulated  and  compared  at 
these  two  angles  to  see  if  any  differences  could  be 
distinguished.  However,  for  a  constant  delay  setting, 
the  two  shape  spectra  were  found  to  be  identical  (see 
fig.  9). 

This  result  is  actually  not  too  surprising.  It  can 
be  shown  that  the  pulse  shapes  formed  in  a  thin  detector 
tend  to  be  smeared  out  because  of  the  finite  electron 
scattering  radius,  and  hence  any  subtle  shape  differences 
which  would  occur  due  to  a  polarization  effect  could 
be  destroyed.  An  example  of  this  effect  will  be  discussed 
later.  What  is  required  for  good  shape  measurements  is 
a  reasonably  thick  planar  detector,  in  which  the  average 
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Delay  a  =  9  ns  a  =  50  ns 

Figure  9  Experimental  Delta  Spectra  For  A  Thin 
Planar  Detector 


photon  scattering  distance  is  much  larger  than  the 
electron  scattering  radius  R0.  However,  since  no  such 
detector  was  available  for  experimentation  in  the  lab, 
further  experimental  results  had  to  be  simulated  using 
Monte-Carlo  computing  techniques. 

The  results  of  these  computer  programs  are  shown 
in  fig.  10.  These  graphs  illustrate  the  one-compton 
$  spectrum  for  a  detector  of  radius  1.73  cm  and  thickness 

O  Q 

3.4  cm,  oriented  at  <*)  «  0  and  90  .  A  0.4  MeV  gamma 
source  with  a  polarization  of  +1  was  assumed,  and  all 
pulse  shapes  were  considered  to  be  sharp.  The  statis¬ 
tical  error  in  these  graphs  is  estimated  to  be  110#. 

It  should  be  noted  that  the  actual  b  spectrum  for  this 
experiment  would  include  other  contributions  from  both 
photoelectric  and  multicompton  events,  and  since  these 
interactions  are  for  the  most  part  polarization  insen¬ 
sitive,  they  will  tend  to  decrease  any  polarization  effect 
inherent  in  the  one-compton  &  spectrum. 

These  graphs  suggest  that  £  is  not  very  sensitive 
to  the  small  polarization  effect  discussed  above.  This 
result  was  already  hinted  at  in  figs.  7  and  8,  which 
demonstrated  that  &  as  measured  by  a  TCP  is  a  very  poor 
measure  of  the  seperation  of  events  occuring  in  a  one- 
compton  scattering  interaction.  Thus  it  is  extremely 
doubtful  that  polarization  could  be  measured  using  such 


a  device. 
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10  Delta  Spectra  For  A  Thick  Planar  Detector 


U  =  0 

mean  =  5*35 


10 


u)  =  9o° 

mean  =  5-32 


10 


3«2  A  Possible  New  Method  of  Shape  Measurement 
Theory: 

One  problem  encountered  by  the  TCF  was  that  its 
output  was  partially  dependant  upon  the  location  of  the 
first  gamma  interaction  in  the  detector.  It  would  be 
of  an  obvious  advantage  to  design  a  pulse  shape  measure¬ 
ment  device  whose  output  depended  only  upon  the  separation 
of  the  two  gamma  interactions  in  a  one-compton  scattering 
event,  and  not  upon  the  individual  location  of  either 
interaction. 

Fig.  11  shows  the  derivative  of  the  pulse  shapes 
shown  earlier.  It  will  be  noticed  that  the  position  of 
the  first  "kink"  in  this  pulse  shape  occures  at  a  time 
t^  d/2  -  ,  where  d/2  is  the  half  thickness  of  the 

detector,  and  Z ^  the  distance  in  the  z  direction  from 
the  nearest  detector  edge  to  the  outermost  gamma  inter¬ 
action.  Similarly  the  second  kink  occures  at  a  time 
t2  d/2  -  Z2,  where  Z2  is  the  distance  from  the  nearest 
detector  edge  to  the  other  gamma  interaction.  Hence 
notice  that 

€  s  t2  -  ti 
aZ2  “  21  ' 

Thus  an  £  spectrum  of  many  one-compton  pulse  shapes 
should  be  polarization  sensitive. 

Recently  a  new  pulse  shape  analyser  has  been 

2 

developed  which  may  be  capable  of  making  £  measurements.) 
However,  several  factors  dictate  what  an  overall  £ 
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t1  t2 


t^j  <x  d/2  -  Z^ 

Figure  11  Derivative  of  a  one-compton  t£  oc.  d/2  -  Z^> 
pulse  shape . 


Figure  12  Pulse  shape  derivative  assuming  a  sharp 
pulse  shape. 
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spectrum  will  actually  look  like,  and  these  should  be 
considered  before  any  experimental  work  is  carried  out. 

They  are: 

Detector  Thickness  It  is  obvious  that  as  the  thickness 
of  a  detector  is  decreased,  so  to  will  its  maximum  £ 
value  and  average  £  value  decrease.  Therefore  a  thick 
detector  will  have  a  higher  polarization  sensitivity 
than  a  thin  detector  (see  table  1). 

Electron  Range  Fig.  4  showed  the  range  R0  of  an  electron 
of  energy  E  MeV  inside  a  germanium  crystal.  This  finite 
range  will  tend  to  smooth  out  the  kinks  in  the  various 
pulse  shapes  produced,  and  can  as  suggested  earlier 
completely  obliterate  them.  For  example  fig.  12  shows 
the  expected  pulse  shape  derivatives,  without  considering 
electron  range,  for  typical  one-compton  gamma  interactions 
in  two  planar  detectors  of  different  thickness.  It  is 
apparent  that,  if  the  electron  range  effect  were  taken 
into  consideration,  the  pulse  shape  derivative  for  the 
thinner  detector  would  be  less  distinct  than  that  of 
the  thicker  detector.  Hence  the  thinner  crystal  is 
again  found  to  provide  £  measurements  which  are  less 
polarization  sensitive. 

Photoelectric  Effect  Photoelectric  events  within  a 
detector  will  provide  £  measurements  which  are  polarization 
insensitive.  Hence  at  low  energies  where  the  photoelectric 
effect  dominates  over  compton  effects,  the  detector  will 
lose  its  polarization  sensitivity. 
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Detector 

One-Compton  Events 

All  Events 

Thickness 

£mean 

£mean 

Q 

£mean 

£mean 

Q 

(cm) 

w=0 

o)=90 

w=0 

w=90 

(mm) 

(mm) 

(mm) 

(mm) 

0.6 

0. 86±0.04 

0.68±0.03 

0. 12±0.04 

0.78±0.04 

0.72±0.04 

0.04±0.04 

3.4 

3. 7±0. 2 

2.7±0. 1 

0. 16±0.04 

2.9±0. 1 

2 . 5±0. 1 

0.07±0.04 

Table  1  Polarization  sensitivity  in  a  planar  Ge(Li) 
detector  at  a  gamma  ray  energy  of  0,847  MeV. 


One-Compton  Events 

All  Events 

Energy 

(MeV) 

£mean 

w=0 

(mm) 

£mean 

uf90 

(mm) 

Q 

£mean 

co=0 

(mm) 

£mean 

co=  90 

(mm) 

Q 

0.400 

3 . 3±0. 2 

2 . 2±0. 1 

0.20±0.05 

3.4±0.2 

3.1 ±0.2 

0.05±0.05 

0.847 

•  3.7±0.2 

2.7±0. 1 

0. 1 6±0 . 04 

2.9±0.2 

2.5±0. 1  ■ 

0.07±0.04 

1.500 

3.2±0.2 

3.0±0.2 

0.03±0.06 

2.3±0.1 

2 . 4±0. 1 

-0.02±0.04  ■ 

Table  2  Polarization  sensitivity  as  a  function  of  gamma 
ray  energy  in  a  planar  detector  of  thickness 
3,4  cm. 


m 

Multicompton  Events  and  Pair  Production  Since  £  measure¬ 
ments  are  dependent  only  upon  the  two  outermost  events 
in  the  crystal,  £  spectra  from  multicompton  events  tend 
also  to  be  polarization  insensitive.  Thus  at  larger 
energies  (>1  MeV)  where  multicompton  and  pair  production 
events  dominate  over  one-compton  events,  the  detector 
will  again  lose  its  polarization  sensitivity. 

Scattering  Angle  -O-  The  average  scattering  angle  -e* is 
dependent  upon  the  energy  of  the  incoming  gamma,  in  that 
as  the  energy  increases,  the  average  angle  &  decreases. 
This  preference  at  higher  energies  towards  forward 
scattering  will  further  decrease  any  polarization  sensi¬ 
tivity  inherent  in  the  £,  measurements. 

These  then  are  the  major  factors  which  dictate  the 
shape  of  a  particular  £  spectrum.  Other  lesser  factors 
could  be  included  (see  fig.  13),  but  the  arguments 
given  above  represent  the  main  limitations  to  polarization 
measurement  by  this  new  method. 

Results : 

The  results  of  various  Monte-Carlo  computing  programs 
are  summarized  in  tables  1  and  2.  The  computer  simulated 
a  detector  of  radius  1.73  cm  and  thickness  either  0.6  cm 
or  3.4  cm,  and  again  assumed  that  all  pulse  shapes  were 
sharp.  It  should  be  noted  here  that  any  device  capable 
of  making  6  measurements  most  likely  would  not  be  able 
to  distinguish  one-compton  events  from  all  other  inter- 
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Figure  13  Identical  Pulse  Shape  Derivatives 

Note  that  an  £  measuring  device  would  not 
be  able  to  distinguish  between  these  two 
different  one-compton  interactions. 


actions.  Therefore  the  right  side  of  table  2  is  in 
fact  the  best  representation  of  what  real  experimental 
results  would  look  like. 


Chapter  IV 


Conclusions 

It  should  now  be  possible  to  compare  the  various 
polarimeters  mentioned  above.  To  begin  with,  of  the  two 
pulse  shape  measurement  devices  discussed,  the  £  measuring 
device  appears  to  be  the  more  polarization  sensitive. 

This  is  logical  since  its  output  is  directly  related  to 
the  separation  along  a  particular  axis  of  the  two  gamma 
interactions  in  a  one-compton  scattering  event,  and  it 
is  this  separation  which  is  directly  related  to  the 
polarization. 

Next,  a  comparison  can  be  made  between  the  polari¬ 
zation  sensitivity  Q  of  the  £  device,  and  the  sensitivities 
of  the  three  other  conventional  polarimeters  (see  table  3)» 
This  table  shows  that,  at  an  energy  of  0.84-7  MeV,  the  £ 
device  has  a  Q  value  of  about  the  same  magnitude  as  the 
thin  planar  detector,  which  is  of  course  considerably 
less  than  the  sensitivities  of  either  the  two-  or  three- 
crystal  polarimeters.  Furthermore,  whereas  the  £  polari- 
meter  suffers  from  the  same  low  energy  difficulties  as 
the  planar  detector,  (an  inability  to  distinguish  between 
photoelectric  and  one-compton  events),  it  also  suffers 
from  a  few  high  energy  (>1  MeV)  problems  which  the  planar 
detector  does  not  encounter.  These  include  its  inability 
to  distinguish  between  one-compton  and  many-compton 
events,  and  the  increasing  electron  range  R0  which  tends 
to  smear  out  pulse  shapes. 
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Type  Of 

Polarimeter 

Pol  ari  zati  on 

Sensi ti vi ty 

3 

Three-Crystal ) 

0 . 44±0 . 02 

Two-Crystal 

0. 35±0.03 

c 

Planar  Detector) 

0. 08±0 . 05 

£  Measurement 

0 . 07±0 . 04 

Devi ce 

Table  3  Polarization  sensitivity  of  various 

Ge(Li)  polarimeters  at  a  gamma  ray  energy 
of  0.847  MeV. 


Further  experimental  study  could  be  carried  out  to 
verify  the  computer  results  given  above.  However,  this 
would  require  the  design  and  construction  of  a  new  pulse 
shape  discriminator,  and  the  purchase  of  a  new  Ge(Li) 
detector  which  would  be  at  least  as  thick  as  the  one 
simulated  above.  In  light  of  the  comparisons  just  made, 
it  is  my  opinion  that  this  extra  work  and  expense  would 
not  be  meritted.  If  in  the  future  it  becomes  necessary 
to  make  gamma-ray  polarization  measurements,  then  it  is 
my  suggestion  that  a  conventional  two-  or  three-crystal 
polarimeter  be  built. 
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Appendix  1 


Polarization  Effects  In  Photoelectric  Absorption 

One  problem  encountered  when  trying  to  develop  a 
single  crystal  gamma  polarimeter  was  that,  at  low 
energies  photoelectric  absorption  predominates  over 
compton  scattering.  It  would  thus  be  very  useful  if 
any  polarization  sensitivity  inherent  in  the  photoelectric 
effect  could  be  used  to  make  gamma  ray  polarization 
measurements. 

Photoelectric  absorption  can  occur  when  a  gamma 
photon  encounters  an  electron  in  the  vicinity  of  an 
atom's  nucleus.  Virtually  all  of  the  energy  of  the 
photon  is  absorbed  by  the  electron,  although  the  presence 
of  the  nucleus  is  essential  in  order  to  conserve  linear 
momentum,  and  the  electron  is  ejected  from  the  atom  in 
some  preferential  direction.  The  differential  cross 
section  for  the  ejection  of  electrons  from  the  atomic 
K-shell  at  low  energies  is  given  by: 

CosV 

din  l<w  + 

where  k  ~  V  ,  oL  =  , 

C  moC2- 

Z  is  the  atomic  number  of  the  absorbing  material, 

R  is  a  constant,  and 

■©•  and  (j>  are  as  previously  defined  in  fig.  1,  in 
which  the  scattered  photon  is  now  replaced  by  the 
recoiling  electron. 

It  is  obvious  from  this  equation  that  at  low  energies 
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electrons  will  be  preferentially  ejected  in  the  -B-  =  90° 
and  (j)  =  0  direction;  that  is  in  the  direction  of  the 
electric  vector  of  the  incident  radiation. 

At  relativistic  energies,  the  differential  cross 
section  must  be  modified  to  become: 

~  c(  A  -V-  6  Cos*<£) 

dsi 

where  A  =  CoS  0*) 

B=  (  I- -  {*/z)(  l  -§Co4«-) 
and  C  =  (.  '  /  oO  ('"&  )  1  f  S(r\  -  ^604©’')^ 

A  simple  interpretation  of  this  equation  similar  to 
that  given  for  the  low  energy  case  is  no  longer  possible. 
However  it  can  be  stated  that  the  angle  for  which 

the  differential  cross  section  intergrated  over  ^  is  a 
maximum  decreases  as  the  energy  of  the  incoming  gamma 
increases.  As  well  the  asymmetry  ratio  A(©*)  defined  by: 

fr(Q')  _  gCer(^>  =  o°)  -  c4c~(<h  ?o°) 
d  c  +ctcr(^=9o°) 

decreases  for  O  =  _  as  the  energy  increases.  These 

two  relations  are  shown  graphically  in  fig.  14. 

In  trying  to  measure  this  polarization  asymmetry, 
it  is  of  course  impossible  to  directly  detect  in  which 
direction  the  electrons  are  ejected.  At  best  it  is 
possible  to  visualize,  in  some  way,  the  "track"  the 
electron  makes  as  it  scatters  and  loses  energy  by 
ionization  inside  a  detector.  (For  instance,  when  using 
a  planar  Ge(Li)  detector  and  a  pulse  shape  analyser 
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similar  to  either  one  of  the  two  described  above,  it 
is  theoretically  possible  to  measure  the  projection  of 
an  electron's  ionization  track  on  a  line  perpendicular 
to  the  plane  of  the  detector.)  Hence  it  is  essential 
that  photoelectrons  not  only  be  ejected  from  their  atoms 
in  some  prefered  direction,  but  also  that  they  continue 
to  travel  in  relatively  straight  lines  as  they  lose 
energy.  Otherwise  the  electrons  will  "forget"  in  which 
direction  they  were  originally  ejected,  and  any  polariza¬ 
tion  sensitivity  inherent  in  their  ejection  will  be  lost. 

A  simplified  theory  of  electron  multiple  scattering 
can  be  used  to  estimate  just  how  straight  the  electron 
tracks  will  be.  This  theory  assumes  that  the  electrons 
only  undergo  deflection  by  Rutherford  scattering  off 
atomic  nuclei,  and  that  these  deflections  are  small. 

This  type  of  scattering  is  described  by  the  cross  section 

^  =  H  Z *  I  ;  Q- «  I 

<Lsi 

where  p  is  the  electron's  momentum.  After  travelling 
a  distance  d,  a  beam  of  electrons  originally  all  moving 
in  the  same  direction  will  have  a  spectrum  of  directions 
approximately  gaussian  in  shape  and  with  a  standard 
deviation: 

/ei>  =  8TTMZleVf  _ i<  dx 

^  J  tpcjiV- 

o 

where  of  course  it  must  be  kept  in  mind  that  pep  will 
decrease  as  the  distance  travelled  increases.  N  here 
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refers  to  the  number  of  atoms  in  the  detector  per  unit 

volume,  while  K  is  given  by: 

K  **  ln(G-  /6L .  ) 
max  iim^ 

where  0^ax  now  represents  the  maximum  scattering  angle 
in  any  one  collision,  and  the  minimum  scattering 

angle.  In  this  appendix  it  will  be  assumed  that 

^nax  ^  an<^  ^min  21  0*^  so  'kkn'k  K  is  a  constant. 

This  assumption  will  not  greatly  effect  the  final  conclu 
sions  to  the  appendix  due  to  the  insensitivity  of  the 
logarithmic  function. 

At  an  energy  of  1  MeV,  about  96$  of  the  energy  loss 
of  an  electron  is  due  to  ionization.  Hence  the  further 
approximation  is  made  that: 

dr  «  JiT  j 

d  X  d  X  /iors 

_  h (  mocV-T^  _  p* 

l  \  I  *  ( i-£l) /  '  J 

where  T  is  the  kinetic  energy  of  the  electron,  and  I 
the  ionization  energy  of  an  atomic  electron.  If  this 
equation  is  combined  with  the  previous  expression  for 
<e*>,  then  the  change  in  with  distance  travelled 

can  be  calculated  for  any  starting  electron  energy. 

Some  examples  of  these  results  are  shown  graphically  in 
fig.  15-  These  results  can  be  compared  with  fig.  16 
which  plots  electron  kinetic  energy  T  as  a  function  of 
distance  travelled.  In  both  figures  it  is  assumed  that 
the  electron  scatters  inside  a  germanium  crystal. 

One  can  clearly  see  from  these  graphs  that  the 
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Figure  15  Average  spreading  of  electrons  as  a  function 
of  distance  travelled  in  a  Ge(Li)  detector. 


Energy  (MeV) 
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Figure  16  Electron  energy  as  a  function  of 

total  distance  travelled  in  a  Ge(Li) 
detector. 
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electrons  forget  in  which  direction  they  were  originally 
ejected  long  before  they  come  to  a  complete  stop  in  the 
detector.  Thus  it  is  quite  doubtful  if  the  photoelectric 
effect  could  be  used  in  a  Ge(Li)  detector  as  the  basis 
of  any  practical  polarization  measurement  techniques. 


Appendix  2 


Monte-Carlo  Compton  Scattering  Simulation  Program 

A  Monte-Carlo  computing  program  was  written  in 
FORTRAN  IV  to  simulate  the  scattering  of  gamma  photons 
inside  a  planar  Ge(Li)  detector.  The  data  produced  by 
this  program  was  then  used  as  input  to  several  other 
programs,  which  in  turn  produced  the  results  shown  in 
fig.  10  and  tables  1  and  2.  A  brief  explanation  of  this 
program,  followed  by  a  listing  of  it,  is  given  below. 

All  numbers  in  parentheses  refer  to  line  numbers  in  the 
listing. 

In  this  simulation,  a  cylindrical  planar  detector 
of  thickness  DLE  and  radius  R02  is  assumed.  The  detector 
dimensions,  measured  in  centimeters,  can  easily  be 
changed  (see  lines  10, 11, 537) »  but  are  not  actual  input 
variables.  The  program  does  require  as  input  the  initial 
gamma  ray  energy  EN  in  MeV,  a  random  number  seed  ISET 
which  will  be  explained  later,  the  total  number  of  gamma 
photons  NN  which  the  program  will  simulate,  and  the  beam 
polarization  angle  XPH1  in  degrees.  XPH1  corresponds  to 
the  detector  orientation  angle  U  shown  in  fig.  2,  where 
all  gamma  rays  are  considered  to  have  a  polarization 
of  +1. 

The  program  is  structured  in  such  a  way  that  a 
gamma  photon  entering  the  planar  detector  edge  on  must 
have  its  first  gamma  interaction  somewhere  inside  the 
detector  crystal  (56-81).  It  then  determines  whether 
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this  gamma  photon  will  be  photo-absorbed  or  compton 
scattered  (86,87),  and  if  photo-absorption  takes  place, 
returns  to  line  52  to  consider  the  next  gamma  ray  entering 
the  crystal.  If  compton  scattering  occures,  however,  it 
calculates  the  scattering  angles  •©•  (94-106) ,  $(109-120), 
and  p  (139-142)  as  defined  in  fig.  1,  and  as  well  the 
energy  of  the  compton  scattered  photon  (150-154).  It 
then  decides  where  the  next  gamma  interaction  will  take 
place  (157 »  148-165),  and  whether  or  not  this  location  is 
inside  the  detector  (169»  170).  If  the  gamma  ray  has  been 
scattered  outside  the  detector  the  program  ignores  this 
photon  and  returns  to  line  52  to  begin  again.  Otherwise 
the  program  returns  to  line  86  to  ask  if  the  gamma  ray 
will  now  be  photo-absorbed  or  compton  scattered  again. 

This  process  is  continued  untill  the  gamma  ray  is 
either  scattered  outside  the  crystal,  or  totally  absorbed 
in  it.  If  the  photon  is  totally  absorbed,  the  program 
keeps  track  of  where  each  gamma  interaction  takes  place 
(192),  and  how  much  energy  is  lost  with  each  interaction 

(195). 

After  the  specified  number  of  incident  gamma  rays 
have  been  studied,  the  program  outputs  (198)  the  total 
number  of  photons  which  escaped  the  detector  (JJ),  were 
photo-absorbed  (KK) ,  were  compton  scattered  only  once  (KL) , 
were  compton  scattered  twice  (KM),  and  were  compton 
scattered  more  than  twice  (KN).  It  also  outputs  the 
total  number  of  counts  in  the  full  energy  peak  (200), 
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which  is  of  course  just  the  sum  of  KK,  KL,  KM,  and  KN. 

A  number  of  subroutines  are  used  throughout  the 
program.  These  are  outlined  below. 

SUBROUTINE  INRAND  uses  a  random  number  seed  ISET  generated 
by  the  user  to  initialize  a  sequence  of  random  numbers. 

The  subroutine  is  thereafter  entered  from  the  main  program 
at  the  statement  ENTRY  RAND  (215) »  which  generates  a 
random  number  RA  between  0  and  1 . 

SUBROUTINE  DETR  determines  the  range  of  the  gamma  ray  in 
the  germanium  crystal. 

SUBROUTINE  THPR  determines  the  probability  of  a  compton 
scatter  occuring  at  an  angle  O-  for  any  angle 
SUBROUTINE  PHPR  determines  the  probability  of  a  compton 
scatter  occuring  at  an  angle  <p  for  a  specific  angle  O-. 
SUBROUTINE  POPR  determines  the  angle  p  associated  with 
a  particular  gamma  ray  scattered  at  given  angles  ©•  and  ^>. 
SUBROUTINE  COORD  converts  from  the  scattering  coordinates 
0,  <p  and  ^  to  the  detector  coordinates. 

SUBROUTINE  PERC  determines  if  a  gamma  photon  of  a  partic¬ 
ular  energy  will  be  photo-absorbed  or  compton  scattered. 
SUBROUTINE  OUT  determines  if  a  particular  gamma  is  inside 
or  outside  the  detector  crystal. 
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142  CALL  POPR(CET( HUM) , SET (NUM ) , CEP (NUM) ,AL,PI,BE) 

143  CBE(NUM)=COS(BE) 

144  SBE ( NUM ) =S I N ( BE ) 

145  C 
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211  DATA  RNORM  /2 147483643  . 

212  KSET  =MOD ( ! $ET/ 5  ) 

213  I  RAND*! START (KSET+1) 

214  RETURN 

215  ENTRY  RAND(RA) 
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281  XX1=XX1+STEP 

282  GO  TO  15 

283  10  XX1=XX1-STEP 

284  15  CONTINUE 
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316  20  X(2,J)=X(2,J)+(A(J,K)*X(1,K)) 

317  X2 -X ( 2 , 1 ) 

318  Y2=X(2,2) 

319  Z2 -X ( 2 / 3  ) 

320  RETURN 
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345  RETURN 

346  10  I  RUT B0 

347  RETURN 

348  END 
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